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Abstract 


The novel human coronavirus known as SARS-CoV-2 poses a serious 
risk to human health. Regretfully, the US Food and Drug Adminis- 
tration (FDA) has approved very few oral antiviral medications for 
the treatment of COVID-19 patients. In the current study, we have 
explored the topological characterization of the orally bio-available 
SARS MP? inhibitor nirmatrelvir which in combination with ritonavir 
(under the brand name Paxlovid) is recently approved for emergency 
use authorization by FDA. Topological indices are a useful tool in 
chemical graph theory to determine the diverse pharmaceutical, biolog- 
ical and physico-chemical properties of a molecule. In this study, we 
ascertain several well-known degree-dependent topological indices for 
the medication nirmatrelvir directly using their common definitions in 
mathematics and alternatively by utilising M-polynomial after deriv- 
ing M-polynomial of nirmatrelvir. In addition, we plot the obtained 
topological indices and the M-polynomial to comprehend the geometric 
behaviour of them. The outcomes can aid in the investigation of the 
physical characteristics of the recently created medications utilised to 
treat COVID-19. 
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1 Introduction 


Caused by SARS-CoV-2 (Severe Acute Respiratory Syndrome CoronaVirus 2), 
the Corona Virus Disease 2019 (COVID-19) is a RNA virus that is reported 
to cause severe pneumonia in serious cases. The catastrophic outbreak of 
this disease has claimed the life of more than 7 million people globally [51]. 

Many strategies have been developed worldwide in unmatched time- 
line to fight this infectious disease. So far, WHO approved vaccines like 
COVAXIN®, COVISHIELD™, Moderna: mRNA-1273, Ad26.COV2.S 
(Janssen COVID-19), BNT162b2 (Tozinameran / Comirnaty) developed 
by Bharat Biotech International Limited, Oxford University and British- 
Swedish company AstraZeneca (manufactured by Serum Institute of India 
Pvt. Limited), Moderna, Johnson & Johnson, Pfizer-BioNTech, respectively, 
are available as prophylactic measures of COVID-19 [3]. Intravenous ther- 
apeutics like remdesivir, sotrovimab, tocilizumab, tixagevimab, etc, were 
found to be effective to cure the infection brought on by the fatal virus [47]. 
The need of oral drugs has now become imperative, as it is believed that 
oral drugs can significantly reduce hospitalization and death. Oral antiviral 
drugs such as Molnuperavir (RNA dependent RNA polymerase inhibitor), 
Paxlovid (nirmatrelvir co-packaged with ritonavir) and AT-527 are currently 
under clinical trials [38, 48]. 

The current study involves the antiviral medication nirmatrelvir (other 
familiar names are PF-07321332, UNU-7R9A5P7H32) which is a covalent 3CL 
protease inhibitor and is developed by Pfizer by modifying Lufotrelvir (PF- 
07304814) which is another intravenous clinical candidate used for Corona 
Virus [39]. The SARS-CoV-2 genome encodes four structural proteins and 
two polyproteins. The polyproteins are cleaved by endopeptidase main 
protease” (MP?°) into short non-structural proteins that are important for 
viral replication. Corona Virus MP" is a cysteine protease having a catalytic 
dyad of Cys145-His41 in the active site. The presence of Gln residue has 
been reported to be present in the substrate of MP’°. As no human cysteine 


5M-pro is also known as 3CL protease. 
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protease is reported to cleave after Gln, MP’° can be used as a potential, 
selective tractable target for designing drugs against Corona virus. According 
to the reports, nirmatrelvir inhibits the action of these 3CL proteases, 
eventually inhibiting the replication of the virus [38]. 

The molecular formula of nirmatrelvir is C23H32F3N5O04. Figures 1(a) 
and 1(b) depict 3-dimensional images of nirmatrelvir with and without hy- 
drogen, respectively, whereas Figure 2 represents the 2-dimensional chemical 
structure of nirmatrelvir [37]. 


(a) A 3-D image of nirmatrelvir (b) A 3-D image of nirmatrelvir 
with hydrogen. without hydrogen. 


Figure 1: 3-Dimensional images of Nirmatrelvir / PF-07321332 / Unii- 
7R9OA5P7H32. 


Figure 2: Two-dimensional representations of the chemical structure of 
nirmatrelvir / PF-07321332 / Unii-7R9A5P7H32. 
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To comprehend a recently introduced drug’s safety profile, effectiveness, 
and tolerability on the human body, a wide number of research experiments 
is acommon practice. Furthermore, these research experiments essentially re- 
quire huge time, extensive laboratory experiments, vast financial expenditure 
and high-tech laboratories, which are exorbitant for economically developing 
nations. Researchers have shown exceptional interrelations between a large 
number of chemical molecules’ experiments and the chemical graphs that 
correspond to them to comprehend their pharmaceutical behaviors, biological 
activities and physico-chemical properties. The above description of the 
drugs can be substantially predicted in the field of chemical graph theory, 
which links mathematical chemistry and graph theory, by avoiding the lab 
experiments [13, 34, 43, 45, 23]. 

As the topological characterization of nirmatrelvir / PF-07321332 / 
Unii-7R9A5P7H32 is not investigated so far, our goal is to investigate the 
different degree-dependent topological indices in order to better understand 
its related properties and pharmacokinetics. Also, it can open a path to 
explore further the correlated drugs employed in treating COVID-19 patients. 


2 Basics and Survey of Relevant Literature on 
Topological Indices and M-polynomial 


Let A = (V,£) denotes a connected, simple, unweighted and undirected 
graph.Its non-empty vertex set and edge set are denoted by V = V(A) and 
E = E(A), respectively. The notation d(w) is used to represent the degree 
of a vertex w € V. Moreover, we use the notation 6 to designate the set 
min{d(w)|w € V(A)}, and A to designate the set max{d(w)|w € V(A)}. 
Usually, a molecular or chemical graph is used to narrate a molecule of a 
chemical compound. Here, the edges and vertices of the chemical graph 
are furnished by the atoms and the chemical bonds that connect them, 
respectively. 

A molecular or chemical graph’s topological index is a real exponent 
which is a notable and efficient tool during data analysis of the chemical 
graph and it can reflect different chemical, biological and physical features 
of that molecular compound. Diverse topological indices are catalogued 
with reference to several topological behaviors of a chemical molecule. In 
fact, quantitative structure-activity relationships (QSAR) and quantitative 
structure-property relationships (QSPR) can actually be found using these 
individually or combined indices. The researchers are keen to deduce the 
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topological behavior of the molecular graph by examining the QSPR/QSAR. 
For detailed discussion, please refer to [13, 34, 43, 45, 23, 31, 32]. 

Usually, these standard topological indices of the chemical graphs are 
categorized into counting related topological indices, spectrum-based topo- 
logical indices, distance-based topological indices, degree-based topological 
indices, and degree-distance-based topological indices. Degree-based topo- 
logical indices are important in the fields of theoretical chemistry, biology, 
and pharmacology. Some of the standard and citable degree-dependent 
topological indices and their usual mathematical definitions are summarized 
in Table 1. 


Table 1: Degree-dependent topological indices with their associated deriva- 
tion formulas via M-polynomial of a graph A. 


Sl. | Topological Proposed : eT Derivation formula 
No.| Index arouad Formula of topological indices f(a,b) over (M(Azi.y)) 
a, | Pitst Zagreb In- | 4 979 M(A)= YD (d(v) +4(w)) a+b (Dz + Dy)(M(A32,4))|e=y=t 
dex [20] vweB(A) 
Second Zagreb 
2 i = Bi 7x, =y= 
b | Index [20] 1972 M,(A) were ae) ab (DzDy)(M(A; 2, y))|e=y=1 
Modified Sec- 
Cc. ond Zagreb In- | 2004 ™M(A)= YO Aaya) 4 (S2Sy)(M (A; 2, y))|2=y=1 
dex [33] vweE(A) : 
General 1975 
d. Randié In-| joao’ Ra(A) = d(v)d(w))* ab) DD?) (M(A; 2, y))|e=y=1 
1998 7 i 
dex (42, 2: 1 : vweE(A) 
nverse Randié | 1975, _ 1 Ya ae) (7 . 
“ ndex [42, 2, 1] | 1998 RRo(A) any (d(v)d(w))™ (ab)* (S55 )(M (Asa, 9) )le=y=a 
Symmetric 5 \) : 
canal mini ,d(w 24 p2 ; : 
fs Division (Deg) | 2010 Z my 12), ee )) ate (DzSy + DySz)(M(A;2,y))le=y=t 
Index[49] vweB(A) nax(d(y). df w)) 
min(d(v), d(w)) 
Harmonic . wy 2 ¢ oy 
99 = J (M(A; = 
5: ndex [18] 1928 H(A) vweB(A) d(v) + d(w) a+b 25, J(M(A; 2%, y))|e=1 
nverse Sum d(v)d(w 
h. | (Indeg) — In- | 2015 ISI(A) = Che ao) at S2JDzDy(M(A;a,y))|2=1 
dex [46] vweE(A) 
P Augmented Za- AZ(A) = 3 ab 3 13, 3 3 () : 
ts greb Index [19] 2010 > d(v)d(w) (g +b— 3) S7Q-2FD,Dy(M(A; ,Y))|a=a 
owEB(A) d(v) + d(w) — 2 


where, Dy(f(x,y)) =% af (x, ¥)) Dz(f(x,y)) =2« (f(2,y)) Sy(f(a,y)) = Nh £(@.8) ap 


S.(f(0,y)) = fe LEMat, Qalf(ey)) =2°F(e,y), @#0, — I(Fl@,y)) = Fla,2). 
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As an alternative for computing the indices of a molecular graph struc- 
ture using their usual mathematical formulas, a number of graphic poly- 
nomials [21] were introduced in the literature. The fundamental concept 
is to determine a closed-form of an algebraic expression of polynomial of 
the working molecular graph and hence deduce the values of the indices by 
employing some particular differentiation or integration operations (or a 
combination of both) over the algebraic expression of the polynomial. 

Some of the well-known polynomials are the matching polynomial [17], 
the M-polynomial [16, 11, 4, 12], the Tutte polynomial [30], the Schultz 
polynomial [22], the Hosoya polynomial [24], the Clar covering polynomial 
(usually known as Zhang-Zhang polynomial) [52], etc. Among them, the 
M-polynomials plays a significant role in finding the degree dependent 
topological indices. The concept of M-polynomial was proposed in 2015 by 
Deutsch and Klavzar [16]. 


Definition 1 (16]) The mathematical expression 


Mi he oy) = SS tiipal A) cPy? 
d<psqsA 


is commonly known as the M-polynomial of a given graph A, where mp q(A) 
denotes the number of edges vw € E(A) having d(v) = p>1 andd(w) =q>1. 


As justified in [15], any degree-based topological index for a graph 
A = (V, E) is symbolized by D(A), and is expressed as 


D(A) = D7 F(dw),d)), 


uve 


where the function f(a, b) is catalogued in Table 1 and in fact, it can also 
be re-written as 
D(A) = > Mp,q(A)F(P, @)- 
PS 
Moreover, Table 1 summarizes all the algebraic relations between the degree- 
dependent indices and their corresponding derivation formulas over the 
M-polynomial of a given graph A [16]. 

To find the M-polynomials and the correlated degree-based topological 
indices for various molecular graph / network structures, extensive research 
investigations have been conducted. Some of the such well-known molecular 
graph / network structures are Coronoid systems [27], Silicon-Carbons [5], 
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Benzenoid systems [28], Molnupiravir [14], polyhex nanotubes [36], Kagome 
lattice [25], convex polytopes [44], different Hex-derived networks [29, 6, 7, 8, 
40, 10, 9, 41], hyaluronic acid-doxorubicin [26], Remdesivir compound [50], 
Chloroquine, Remdesivir, Hydroxychloroquine and Theaflavin [35], etc. 


3 Methodology 


In our current study, we have assessed the molecular structure of the oral 
antiviral medication nirmatrelvir using standard degree-dependent topo- 
logical indices. Figure 2 shows the drug’s two-dimensional chemical graph 
structure. we consider the hydrogen suppressed chemical graph of the 
drug nirmatrelvir, because the vertices corresponding to the hydrogen 
atoms have no effect on the graph isomorphism. See Figure 3. Let us 
denote the graph shown in Figure 3 by the symbol N = (V(N), E(N)) = 
(V (PF-07321332), E( PF-07321332)). 


Figure 3: Graphical structure representation of the hydrogen suppressed 
nirmatrelvir / PF-07321332 / Unii-7R9A5P7H32 with edges and vertices. 


To start with, we compute the considered degree-dependent topological 
indices of the graph N by their respective mathematical formulas stated in 
Table 1, in Section 4. On the other hand, in Section 5, we obtained a closed- 
form of an algebraic expression of M-polynomial by partitioning the edge 
set of the graph N and therefore derived the topological indices by applying 
the respective mathematical operations (which reported in Table 1) over 
the M-polynomial. Additionally, in Section 6, we illustrated the obtained 
M-polynomial and the topological indices (in increasing order) by using 
Maple-13 software by plotting the surfaces, to view their geometrical nature. 
In Section 7, we finally come to a conclusion. 
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4 Definition-based Computation of the Degree- 
dependent Topological Indices of Nirmatrelvir 


We have denoted N = (V(N), E(N)) = (V(PF-07321332), E( PF-07321332)) 
to represent the molecular graph of nirmatrelvir / PF-07321332 which is 
depicted in Figure 3, It follows from Figure 3 that the cardinality of the 
vertex set V(NV) is given by |V(N)| = 35 and the cardinality of the edge 
set E(N) is given by |E(.N)| = 37. In addition, the minimum degree and 
maximum degree of a vertex of N are given by 6 = 1 and A = 4, respectively. 
Let us now partition V(V) of N into four disjoint sets, on the basis of degree 
of the vertices, as 


Vi(N) = {2 € V(N) : d(z) = 1}, 
Vo(N) = {z € V(N) : d(z) = 2}, 
V3(N) = {2 € V(N) : d(z) = 3}, 
and V4(N) = {z € V(N) : d(z) = 4}, 


and their cardinalities are given by |Vi(N)| = 13, |Vo(NV)| = 8, |V3(V)| = 11 
and |V4(.V)| = 3. Moreover, let us partition E(N) of N into seven disjoint 
distinct sets based on the degree of the end vertices of all the edges and 
distinguish them as 


Egrsy = {€ = v102 € E(N) : d(vi1) =r, d(v2) = s}, 


where {r,s} = {1,2}, {1,3}, {1, 4}, {2, 2}, {2,3}, {3,3} and {3,4}, and the 
cardinalities of them are given by |Fy9;| = 1, |Fa33| = 4, |r] = 8, 
|Ey2,23| = 2, |Eye,33| = 11, |Eys,33| = 7 and |Ey3,4)| = 4. 

Let us now compute the degree-dependent topological indices of our 
considered molecular graph of nirmatrelvir, by their respective mathematical 
formulas listed in Table 1. 


Theorem 2 Let N = (V(N), E(N)) = (V(PF-07321332), E( PF-07321332)) 
be the molecular graph of nirmatrelvir. Then 


(i) M,(N) = 192, 
(ii) Mo(N) = 231, 
(iti) ™ Ma(N) = 322, 
(iv) Roa (N) = 2%+4x 3%4+8 x 4%42x 4% +11 x 694+ 7x 9% +44 x 12%, 
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(vy) RRAN)= t+ A+ B+4+bt+ +7 


(vi) SDD(N) = 100, 


(vii) H(N) = 238 
(viti) ISILN) = 33h, 


. _ 60370561 
(ix) AZ(N) = 216000 ° 


Proof: 
their values for graph N are calculated underneath. 
(i) First Zagreb Index: 
M(N)= > (dv) +d(w)) 
vweE(N) 
= 1(14+ 2) 4+ 4(14+3) + 8(1 +4) + 2(24 2) +11(2 +3) 


+ 7(3 + 3) +4(3 + 4) 
= 192. 


(ii) Second Zagreb Index: 


vweE(N) 
= 1(1 x 2)+ 4(1 x 3) + 8(1 x 4) + 2(2 x 2) +.11(2 x 3) 


+ 7(3 x 3) + 4(8 x 4) 
= 231. 


(iii) Modified Second Zagreb Index: 


1 
™ M>(N) = aoe 
: "2n d(v)d(w) 


Wt r Be sy a A a, i 4 
RS ao i ee Oe es ee oe 
131 


= — 7.2 : 
18 T.2777777 


121 


By using the mathematical definition of the topological indices, 
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(iv) General Randié Index: 
Ro(N)= SY) (d(v)d(w))* 
vweE(N) 
= 1(1 x 2)° + 4(1 x 3)* + 8(1 x 4)® + 2(2 x 2)% 4+ 11(2 x 3) 
+ 7(3 x 3)* +4(3 x 4)° 
=2%+4x 3%4+8x« 4%*+2x4%411 x 6%°4+7x 9% +4 x 12%. 


(v) Inverse Randié Index: 


1 
RRAN)= Game 


vweE(N) 
1 4 8 2 11 
~ axa t xa Axa * @xde* @xaye 
ae 
' (3x 3)% ' (3x 4) 


1 4 8 2 11 7 4 


(vi) Symmetric Division (Deg) Index: 


- min(d(v), d(w)) ; max(d(v), d(w)) 
Se aan ae) : vin(ateytayy 


_,f1,2 alas alas Fe 0 
ain eh aa Bae ae | eae a oS 
D2 28 aes Bey al 

oe repay pe! 


“isa tsa} rea} aaa] 
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aed iT a oT : 1 4 e 
28 Bey amen eee 


_ 1548 


~ 105 


516 
= —— & 14.7428571. 
35 


(viii) Inverse Sum (Indeg) Index: 


d(v)d(w) 
ISI(N) = BoC lc 
= d(v) + d(w) 


ay 1x: 44 1x3 Ry ema Ay kes, 
Labo ase tea S| Garo 
2x3 3x3 3x4 
+ 11 + 7 4< —_— 
es} ioah+ eo} 


8951 
= 42.6238095. 
(ix) Augmented Zagreb Index: 
dv)aw) 
AZN)= So { 
ea d(v) + d(w) — 2 


ng lM Veg J ES one ice an le 
Basal (key mee eae in ee eae Ppa 9 
aan ea 2 : 
Sp) +11 ee 
942-2 ia: eae 
ee al las a4 
SU 2: fd 
Hees} taf} 


60370561 
= ——— #8 279.4 : 
216000 79.4933379 


Remark 3 From Theorem 2, we see that Ra(N) % 90.14321042 and RR.(N) + 


15.99527292 for a = 3. 


5 M-polynomial of the Graph of Nirmatrelvir / 
PF-07321332 


We now determining a algebraic expression of M-polynomial for the graph 
structure of nirmatrelvir and subsequently extort the numerical values of the 
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related topological indices (degree-dependent) by using the listed (Table 1) 
derivation formulas. Therefore, we will validate the values of the indices 
calculated in the previous section. 


Theorem 4 Let N = (V(N), E(N)) = (V(PF-07321332), E( PF-07321332)) 
be the molecular graph of nirmatrelvir. Then the expression of M-polynomial 
of the molecular graph is given by 


M(N;2,y) = cy? + 4ary? + 8xy* + Qr7y? + 1127y? + 7x3 y? + 43 y4. 


Proof: We have analysed in the beginning of the previous section that 
the edge set E(.NV) has the cardinality 37 and the edge set can be divided 
into seven discrete subsets based on the degree of the end vertices as 


Eg sy = {e = v102 € E(N): d(v1) = 7, d(v2) = s}, 


where {r,s} = {1,2},{1,3} {1,4}, {2,2} {2,3}, {3,3} and {3,4}, having 
frequencies given by |F49}| = 1, |Fr33] = 4, |Fa| = 8 |Ey223| = 2, 
|Eyo.33| = 11, |Fy3,33| = 7, and |£y3.44| = 4. Therefore, for our considered 
graph N, using the formal definition of the M-polynomial, we have 


M(N;2,y) = De Mp g@?y?, where p,q € {1,2,3,4} andéd=1,A=4 
d<p<qsA 
=mpry’ +m 3ry? + myasry! + mo227y" 
+ mo3r7y° a mM3,30°Y° - m340°y" 


= |Ex 2} cy” + [E4331 ay? + |Eq yl cy + |E'¢2,93| ay? 


+ |Ep33| ay? + |E¢3,3} | ayer |E¢3.44| ay? 
= ry” + doy? + 8ry* + 2a7y? + 1127 y3 + Tr3y? + 4a3y?, 
The following theorems extract the values the previously discussed topo- 
logical indices of the graph N from our obtained M-polynomial M(N;z, y). 
The reader gets to see some of the proofs for the following lemmas. 


Lemma 5 Let N be the graph of nirmatrelvir, then the differential operators 
Dz =a, Dy =’, and D;Dy over the M-polynomial M(.N; x,y) are given by 


Dz (M(N; 2, y))=xy? + 4ay?+ 8aryt+ 4a7y?4+ 2977 y3 + 21a? y?+ 1207 y4, 


Dy (M(N; 2, y)) =2cy? + 12ay? +320y4+407y?+3327y3 + 2123 y3 +1627 y4, 
¥y 
DzDy(M(N; 2, y)) =2ry*? + 1l2ry? +32zry* +8a7y? + 66a7y? +63a%y? +4827 y4, 
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Theorem 6 Let N be the graph of nirmatrelvir, then the first and second 
Zagreb indices of the graph N are given by Mi(N) = 192 and M2(N) = 231, 
respectively. 


Proof: Using Lemma 5 and the main derivation formulas mentioned in 
Table 1, we have 


My(N) = (De + Dy)(M(N32,y))|e-yo1 = 72 + 120 = 192. 


M2(N) = DrDy(MN; 2, y))|e=y=1 = 231. 


Lemma 7 Let N be the graph of nirmatrelvir, then the integral operator 
SzSy over the M-polynomial M(N; x,y) is given by 


i 4 1 fl 7 1 
Sy Sy(M(N;2,y)) = sey" + ,2y?+20y!+ 50 + ay’ + ow yh + 2a y!. 


Proof: Let us first impose the integral operator S, over the M-polynomial 
M(N;«x,y) and then S,, on the obtained result. 


M(N; x,t) 
f 


 $,(M(N;2,y)) = iH at 


1 4 ‘i 7 

= Lay? + Lay? + Qayt + 2?y? + Kay? + Lo3y3 + 23 y!, 
2 3 5 3 

and 

SxSy(M(N; x, y)) = Sx(Sy(M(N; x, y))) 


1 4 1 11 7 1 
— sty + gty' t+ 2ay'+ ins aout ities ay. 


Theorem 8 Let N be the graph of nirmatrelvir, then the modified second 


Zagreb index of the graph N is given by ™M2(N) = z. 


Proof: | Using Lemma 7 and the main derivation formula mentioned in 
Table 1, we have 


eee = T.2777777. 


™M2(N) = SxSy(M(N; 2, y))|e=y=1 -_ 18 
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Lemma 9 Let N be the graph of nirmatrelvir, then the differential operator 
Dy Dy and the integral operator SYS7 over the M-polynomial M(N;2,y) are 
given by 


Dz Dy (MN; x, y)) 
= gy? +4 x 3%xy? +8 x 4%ayt +2 x 2°2% ey? +11 x 2°3%r7y3 
+7 x 3%3%r3y3 + 4 x 374% 23 y4, 
SxSy(M(N; x, y)) 
1 9 4 gy 8a, 2p oy os 


= 5aty Tay d wy t saqat Ut 


4 
30 4a ny, Jaga 


ay’. 


3,3 
Theorem 10 Let N be the graph of nirmatrelvir, then the general Randié 


index Ra(N) and inverse Randié index RRq(N) of the graph N are given by 


Ra(N) = 2%°+4 x 3%4+ 8 x 4%+2 x 4%4 11 x 674-7 x 9°+ 4 x 12%, 


‘ee eee: omen eee 4 
and RR(N) = 35 3a + gat gat gat gat JQ" 


Proof: Using Lemma 9 and the main derivation formulas mentioned in 
Table 1, we have 


Ra(N) = DE Dy (MN; #, y))le=y=1 
= 2°+4 x 3°+8 x 4%+2 x 4%4+ 11 x 6°+7 x 9°44 x 12%, 


and RR(N) = S7S7(M(N; 2, y))|e=y=1 
eee ee ee ee eee 


Lemma 11 Let N be the graph of nirmatrelvir, then the operators Dz Sy 
and DySz over the M-polynomial M(N; x,y) are given by 


1 4 22 
D,Sy(M(N; 2, y)) = sey + ary + Qry* + 2a7y?4 3 ey + Tx? y>+ 30% y4, 


and 


33 16 
DyS;,(M(N;2,y)) = Qry” + 12xy?+32ry*+2x7y"4 5 wey’ + Tay + ays. 


Theorem 12 Let N be the graph of nirmatrelvir, then the symmetric divi- 
sion (deg) index of the graph N is given by SDD(N) = 100. 
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Proof: Using Lemma 11 and the main derivation formula mentioned in 
Table 1, we have 


SDD(N) = (DrSy + DySz)(MUN; 2, y))|r=y=1 = 69 + 31 = 100. 


Lemma 13 Let N be the graph of nirmatrelvir, then the operator S;,J over 
the M-polynomial M(N;2,y) is given by 


1 3 19 4 
SaJ (MN; 2,y)) = rc + 5t + 5 2° a’. 


ve 


pt gh 
6 


Theorem 14 Let N be the graph of nirmatrelvir, then the harmonic index 
of the graph N is given by H(N) = ae 
Proof: Using Lemma 13 and the main derivation formula mentioned in 
Table 1, the harmonic index of the graph N, derived from the M(N;z, y), is 
calculated as 


HUN) = 28,J(M(V;2,y))lona =2(5 +}~+—+4-4 


i) 
Ol 
o>) 
NI 


1548 516 
= — = — & 14.7428571. 
105 39 


Lemma 15 Let N be the graph of nirmatrelvir, then the operator S;JD,Dy 
over the M-polynomial M(.N; x,y) is given by 
98 21 48 

z° 4 Ba x. 


T oe 


5 2 7 


2 
SgJD,D,(MN 2,9) = rl + 504 + 


Proof: We know that J(f(z,y)) = f(x,x). Using the computed term 
D,D,(M(N;2,y)) in Lemma 5, we have 


SJ DzDy(M(N; 2, y)) 
= S,J(2ry? + 12axy* + 32xy* + 8x7y? + 66x7y? + 63x%y° + 48a%y*) 
= $,(22° + 20a* + 98x° + 632° + 482") 
2 98 21 48 
=30 + 5a 4 ao ie: 7a. 
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Theorem 16 Let N be the graph of nirmatrelvir, then the inverse sum 


. a: _ 8951 
(indeg) index of the graph N is given by ISIN) = y5- 


Proof: Using Lemma 15 and the main derivation formula mentioned in 
Table 1, the inverse sum (indeg) index of the graph N, derived from the 
M(N;2,y), is calculated as 


ISI(N) = S;JDzDy(M(N; 2, y))|2=1 


2,5, 98, 21, 48 _ 8951 
SB eS ae So TG 


& 42.6238095. 


Lemma 17 Let N be the graph of nirmatrelvir, then the operator S2Q-3/ DED; 
over the M-polynomial M(.N; x,y) is given by 
S:Q-2JD;Dy(M(N; 2,4) 

BAe es OES 
a 33 xu + 38 x 
7x3? x3? 4) 4x3? x4? 5 
T E 4 ba 53 oe 
Proof: We know that J(f(zx,y)) = f(x,xz) and Qa(f(z,y)) = 2° f(z, y), 
a #0. Let us first apply the differential operator Dap; on our obtained 
M-polynomial M(N;2,y) and then impose the remaining operator S3Q_9J 
on the attained result: 

S7:Q-2JD; Dy (M(N; 2,4) 

= S8Q-2J (2? zy? +4x 33ry3 +8 x Bary* + 2 x 23 x 2377y? 

+11 x 2° x 382?y3 +7 x 3° x 382%y? +4 x 33 x 43r%y*) 
= S3Q_9(28a* +4 x 382t + 8x 4a° +2 x 23 x Bat + 11 x 23 x 332° 
+7 x 3° x 382° +4 x 3° x 432") 

= So ( eax srs x a +22" Re 11 x2? 6 328" 

+7x 38 x 3824 + 4 x 33 x 482°) 
BST 5 VORA 2 DID? 

ae ay 38 Ya 38 x 

PROS ay CRG RA ce 


ABR ae gees ae 
aa 


4 x 33 
=— Bert 8 x 


PLO ea 
x 


_ 93 
= 2°a+ 38 
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Theorem 18 Let N be the graph of nirmatrelvir, then the augmented Zagreb 


index of the graph N is given by AZ(N) = Oe 


Proof: Using Lemma 17 and the main derivation formula mentioned 
in Table 1, the augmented Zagreb index of the graph N, derived from the 
M(N;2,y), is calculated as 


AZ(N) = S3Q-2JD3D3(M(N;2,y))|2=1 


3)° Ay? 9)? i? 
= 22 442° = 23690? 4.11 OP be Fe ne game 
va{Sh eft paxaenxvsr{ th pal 
60370561 
= ey 970,4039970, 
216000 


Therefore, we see that with the help of the derivation formulas (as listed 
in Table 1) over our proposed M-polynomial M(N; x,y) of the graph N, it 
is an alternative, quick and compact way of finding the standard degree- 
dependent topological indices. 


6 Geometrical Interpretation of the M-polynomial 
of Nirmatrelvir 


To view the geometrical characteristics of the M-polynomial of nirmatrelvir, 
we draw the 3-Dimensional surfaces of the M-polynomial of nirmatrelvir by 
the Maple-20 in numerous domains, such as, —0.5 < x,y < 0.5 in Figure 4(a), 
—1<z,y <1 in Figure 4(b), -2 < 2z,y < 2 in Figure 4(c), -3< z,y <3 in 
Figure 4(d), —4< x,y < 4 in Figure 4(e) and —5 < x,y < 5 in Figure 4(f). 
These surfaces seem to show similar behaviors but in reality, they have 
different gradients in different regions. Note that, the parameters x and y 
allow us to regulate the M-polynomial value. 


(a) 3-D surface of the —_(b) 3-D surface of the — (c) 3-D surface of the 
M-polynomial of M-polynomial of M-polynomial of 
nirmatrelvir in nirmatrelvir in nirmatrelvir in 
—0.5 <a2,y < 0.5. -1l<ay<1. —2<a24,y <2. 


130 Shibsankar Das and Modhuleena Mandal 


2 > 1 2 4 2 2 
a 2 Papa 2 44 2 


(d) 3-D surface of the —_(e) 3-D surface of the — (f) 3-D surface of the 


M-polynomial of M-polynomial of M-polynomial of 
nirmatrelvir in nirmatrelvir in nirmatrelvir in 
=4 <7. se o. —A<ay<4. =o & ty A 8: 


Figure 4: The 3-Dimensional geometrical illustration of the M-polynomial 
of nirmatrelvir in numerous regions of x and y. 


In order to observe the comparative values, the maximum, minimum 
and intermediate values of the degree-dependent topological indices of nir- 
matrelvir, we plot them, in increasing order, in Figure 5. 


300 


150 


100 90.14321042 


Values of Topological Indices 


14.7428571 15.99527292, 


7.277177 


Modified Harmonic Index Inverse Randié Inverse Sum General Randié. Symmetric “First Zagreb Second Zagreb Augmented 
Second Zagreb Index (a=1/2) _(indeg) Index Index (a=1/2) _ Division (deg) Index Index Zagreb Index 
Index Index 


Different Topological Indices 


Figure 5: Visualizing the values (in increasing order) of the standard degree- 
based topological indices of nirmatrelvir. 


For the molecular graph N = (V(N), E(N)) = (V(PF-07321332), 
E( PF-07321332)) of the drug nirmatrelvir, the degree-dependent topological 
indices are listed in ascending order as follows: modified second Zagreb 
index, harmonic index, inverse Randi¢ index (@ = 1/2) , inverse sum (indeg) 
index, general Randi¢ index (a = 1/2), symmetric division (deg) index, first 
Zagreb index, second Zagreb index, augmented Zagreb index. 
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7 Conclusions 


In the current study, we have considered a bio-available oral antiviral drug 
nirmatrelvir (also familiar as PF-07321332, UNII-7R9A5P7H32) which is 
designed and developed by the scientists of Pfizer. Nirmatrelvir co-packaged 
with ritonavir tablet is available in the brand name of Paxlovid which is 
recently granted for emergency use authorization by U.S. FDA to fight 
against the Corona Virus. Here, we have analyzed some topological at- 
tributes of nirmatrelvir which will aid in perceiving the pharmacokinetics 
and other characteristics of the drug. Some well-known degree-dependent 
topological indices of nirmatrelvir were calculated by their corresponding 
mathematical formulas. We also recalculated and validated the indices from 
the M-polynomial of nirmatrelvir, which is a very swift and precise way to 
compute the degree-dependent topological indices. The obtained numerical 
values of the indices and the M-polynomial were assessed graphically to 
understand their geometrical behavior. This study can serve as new av- 
enue for the investigation of related drugs employed in the treatment of 
COVID-19 patients. 
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